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(-H ' An analytical model founded on geometric and potential energy principles for kink band de- 

formation in laminated composite struts is presented. It is adapted from an earlier successful 
study for confined layered structures which was formulated to model kink band formation in 
the folding of geological layers. The principal aim is to explore the underlying mechanisms 
governing the kinking response of fiat, laminated components comprising unidirectional com- 
posite laminae. A pilot parametric study suggests that the key features of the mechanical 
response are captured well and that quantitative comparisons with experiments presented in 

■ the literature are highly encouraging. 
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1 Introduction 

Kink banding is a phenomenon seen across many scales. It is a potential failure mode for any 
' layered, laminated or fibrous material, held together by external pressure or some form of internal 

$H . matrix, and subjected to compression parallel to the layers. Many examples can be found in 

the literature concerning the deformation of geological strata [H [HIS] , wood and fibre composites 
[H [3 ini IZl [HI ini EH]) and internally in wire and fibre ropes [TTl [T2]. There have been many 
attempts to reproduce kink banding theoretically, from early mechanical models |13l 114] , to more 
sophisticated formulations derived from both continuum mechanics |15| . finite elasticity theory 
[16| and numerical perspectives for more complex loading arrangements jl7j. 

There has been much relevant work on composite materials with significant problems being 
encountered as outlined thus. First, although two dimensional models are commonly employed 
[181 119] ■ modelling into the third dimension adds a significant extra component. It inevitably 
involves a smeared approach in the modelling of material properties since there is a mix of laminae 
and the matrix with the possibility of voids. Secondly, failure is likely to be governed by nonlinear 
material effects in shearing the matrix material |20j . and this is considerably less easy to measure 
or control than the combination of overburden pressure and friction considered in work on kink 
banding during geological folding [211 HH [23] . 

In the current paper, a pilot study is presented where the discrete model formulated for kink 
banding in geological layers is adapted such that it can be applied to unidirectional laminated 
composite struts that are compressed in a direction parallel to the laminae. This is achieved by 
releasing the assumption that voids are wholly penalized since, in the current case, no overburden 
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pressure actively compresses the layers in the transverse direction. Therefore, the rotation of the 
laminae during the formation of the kink band causes a dilation, which is resisted by transverse 
tensile forces generated within the interlamina region. The coincident shearing of this region also 
generates an additional resisting force, but, as mentioned above, this can be subject to nonlinearity, 
in particular a reduced stiffness that may be either positive (hardening) or negative (softening 
perhaps leading to fracture), which is currently formulated with a piecewise linear constitutive 
law. Work done from dilation and shearing are evaluated; additional features from the original 
model: strain energies from bending and direct compression, and the work done from the external 
load can be incorporated without significant alterations. An advantage of the presented model is 
that the resulting equilibrium equations can be written and solved entirely in an analytical form 
without having to resort to complex continuum models or numerical solvers. 

The primary aim of the current work is to lay the foundations for future research. The geomet- 
ric approach has yielded excellent comparisons with experiments for the model for kink banding 
in geological layers; the same is true currently with the present model being compared favourably 
to previously published experiments [5]. Moreover, the relative importance of the parameters gov- 
erning the mechanical response are also identified in the current study. From this, conclusions are 
drawn about the possible further studies that would extend the current model to give meaningful 
comparisons with the actual structural response for a variety of practically significant scenarios. 

2 Review of model for geological layers 

A discrete formulation comprising springs, rigid links and Coulomb friction has been devised to 
model kink band deformation in geological layers that are held together by an overburden pressure 
[22] . It was formulated using energy principles and key parts of the model are shown in Figures [1] 
and[5] It has been compared very favourably with simple laboratory experiments on layers of paper 




Figure 1: Basic configuration of the discrete model for kink banding in n geological layers. The 
tectonic load on each layer is P with a horizontal reaction force H, the axial stiffness of each layer 
is fc, the overburden pressure on the layers is q and the kink band orientation angle is l3. 

that were compressed transversely and then increasingly compressed axially to trigger the kink 
band formation process. The testing rig used in that study is shown schematically in Figure [Sja) 
and a typical test photograph is shown in Figure ^h) . Assuming that the layers were transversely 
compressible, the kink band orientation angle /3 was predicted theoretically for the first time; it 
being related purely to the initially applied transverse strain derived from the overburden pressure 
q. Figure 2] shows the characteristic sequence of deformation with (a) showing the undeformed 
state with applied overburden pressure and the transverse pre-compression defining /3; (b) showing 
the point where the interlayer friction is released when the internal transverse strain within the 
kink band is instantaneously zero and the band forming very quickly in the direction of /3. It was 
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Figure 2: Two internal layers of the geological model. The kink band width is 6, the normal contact 
force between layers is N, the coefficient of friction is /i, the stiffness of individual rotational springs 
modelling bending is c, the elastic stiffness of surrounding medium per unit layer is fc/ and the 
kink band angle is a. A key assumption is the transverse compressibility of the layers. 

later demonstrated that beyond the condition shown in Figure SJc) , where all the layers have the 
same thickness, whether internal or external to the kink band, lock-up begins to occur as shown in 
FigureSJd). This marked the point where new kink bands formed and these could also be predicted 
by this approach after some modifications were made to the model |23) . An example comparison 
between the theory and an experiment from that study is presented in Figure [5] Moreover, this 
model has also been demonstrated to be suitable for modelling internal kink band formation in 
individual composite fibres under bending that are common in fibre ropes |12) . 

3 Pilot model for laminated composite struts 

As discussed above, the system studied in [22 had layers that were bound together by the mecha- 
nisms of overburden pressure and interlayer friction. The deformation was in fact only admissible 
geometrically if the layers were transversely compressible; the relationship between the kink band 
angle a, which could vary, and the orientation angle (3, which was fixed, being such that interlayer 
gaps, or voids, were not created. For a laminated strut, most experimental evidence from the 
literature also suggests that the kink band orientation angle /3 is basically fixed for each laminate 
configuration 5, 7_. It is noted, however, in a recent study on laminates under combined compres- 
sion and shear that this angle can change as the kink propagates but the angle reaches a limit [24] ; 
in the current work, /3 is taken as a constant equivalent to this limiting value from the beginning 
of the kink band deformation process, which is a simplifying assumption. 

Kink band deformation in laminates involves different mechanisms that incorporate the inter- 
lamina region comprising the laminae and the matrix that binds the component together. Since 
the matrix is itself deformable and that there is no overburden pressure to close any voids, the 
model needs significant modifications to account for the different characteristics of the laminated 
strut. It is worth noting that the assumption for the lay-up sequence of the composite in the 
present case is such that no twisting is generated from the applied compression. Figure [5] shows 
the adapted 2-layer model which omits the following features that are not relevant in the current 
case: the foundation stiffness and the overburden pressure, i.e. q = kf = 0. The kink band 
formation is thus intrinsically linked to the deformation of the interlamina region within the strut. 

Shearing within the interlamina region is the analogous process to sliding between the layers 
in the model for geological folding; the latter being modelled in the energy formulation as a work 
done overcoming the friction force. A piecewise linear model is used to simulate the force versus 
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(a) (b) 



Figure 3: (a) Schematic of the experimental rig used for testing the geological model, (b) A typical 
deformation profile in a physical experiment showing a sequence of kink bands with corresponding 
orientation angles /3. 



displacement relationship in terms of the shear resistance (see Figure [7|) , where fracture modes 
that are relevant for a linear-softening response, see Figure [TJa) , are defined as in Figure [H 

Tensile expansion, or dilation, of the interlamina region is modelled, however, with a purely 
linear elastic constitutive law. In the model described in ij2l it was argued that when the interlayer 
contact force was released that not only would the friction be released but also that the overburden 
pressure would inhibit the formation of subsequent voids within the layered structure. Since in the 
current case there is no overburden or lateral pressure as such, potential dilation of the interlamina 
region needs to be included. As in the previous model, however, the lamina deformation is assumed 
to lock-up and potentially trigger a new band forming when a > 2/3; transverse compression in 
adjacent laminae would then be occurring and stiffening the response significantly. Hence, it is 
reasoned therefore that it would be energetically advantageous for the mechanical system to form 
a new kink band rather than continue to deform the current one [231 ■ 



3.1 Potential energy formulation 
3.1.1 Interlamina dilation 

The resistance to interlamina dilation while the kink bands deform is modelled with a linear 
constitutive law; the dilation resisting force -Fi relating to the dilation displacement ^i, thus: 



Fi{a) = CiSi, Siia) = t 



cos(a — j3) 
cos/3 



- 1 



(1) 



with Ci being the transverse stiffness of the laminate, related to the transverse Young's modulus, 
and t being the thickness of an individual lamina. Since the area over which the interlamina region 
dilates depends directly on the kink band width b, the stiffness Ci can be expressed as: 



Ci = bdki, 



(2) 



where ki is the transverse stiffness per unit area of the laminate and d is the breadth of the strut. 
However, with the lamina assumed to be transversely incompressible in the current model and 
the dilation displacement being assigned purely to the softer interlamina matrix material, a clear 
departure from the geological model, the current lamina thickness is thus t rather than t cos(a — /3). 
This is shown in Figure [5] and is detailed in the highlighted area of that diagram. The relationship 
in equation ([1]) for Si is thus obtained from taking the length AB from FigurelHl where Si = AB — t. 
Hence, there is a transverse tensile strain developed as the gap between the laminae grows as a 
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(c) a = 2/? 



(d) a > 2/3 



Figure 4: Sequence of kink band deformation in the geological folding model, (a) initial state with 
/3 defined from applied q; (b) instantaneous release of contact and hence friction within the kink 
band when a = /?; (c) layers inside and outside of the kink band all have equal thickness when 
a = 2/3; (d) lock-up occurs when a > 2/3 and a new band would form. 



increases from zero to /3. The gap subsequently begins to reduce; when a — 2/3 the gap returns to 
zero, marking the commencement of lock-up. 

The work done in the dilation process is therefore given by the expression: 



Ud 



Fi{a')d\t 



cos(a' — (3) 
cos/3 



kibdi"^ 



cos(a — /3) 
cos (3 



(3) 



It is assumed that the interlamina region would not be damaged in the process of dilation and that 
the only nonlinearity in the constitutive law would be under shear. This is because the dilation 
displacement is relatively smaller than the shearing displacement that is discussed next; this has 
the additional advantage of maintaining model simplicity such that any mixed mode fracture 
considerations can be left for future work. 



3.1.2 Interlamina shearing 

Interlamina shearing or the laminae sliding relative to one another is modelled with a piecewise 
linear constitutive law with the force resisting shear Fu relating to the shearing displacement Su, 
thus: 

Fii(a) = Cii^ii, 6u{a) = [sin(a - /3) + sin/3] , (4) 

cos p 

with Cii being the shearing stiffness of the combination of the matrix and laminae sliding relative 
to one another. The relationship for Su in terms of a and /3 in equation ([4]) is given by examining 
the length BC in Figure |6l However, since the band is basically assumed to form instantaneously 
before any rotation occurs, it is implied that Fii(O) = Su{0) = 0. However, since (3^0, the 
expression = BC +t tan (3 is obtained, where the force and displacement conditions are satisfied. 

When the shearing displacement reaches the initial proportionality limit, i.e. when Su — Sq 
(see Figure [7]), the relationship between Fu and 6u changes to: 

Fu = CuSc ( , (5) 
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Figure 5: Graphs showing the total axial load P (left) and the kink band width b (right) both 
versus the total end-displacement A, taken from [23]. The load versus displacement graph shows 
the lockup of the current kink band as a circle on the theoretical curve. 



where is the shearing displacement when the corresponding resistance force reduces to zero. 
Now, if Su{ac) = Sc and (5ii(q!m) = Sm, the expressions for the resisting force can be written thus: 



Ft, 



Cut [sin(a — /3) + sin /3] / cos /3 for a = [0, ac] 



Cut 

cos/3 



sin(Q; — /3) — sin(Q!M — /?) 



sin(Q;c — 13) — sin(aM — /3) 



[sin(ac — /3) + sin j3] for a > ac 



and a = [ac, "m] if sgn(5c) = sgn(5M) > 0, 
for a ^ au and sgn((5c) = sgn((5M) > 0. 

(6) 

Moreover, since the shear area of contact depends on the kink band width b, the stiffness Cu can 
be expressed as: 



Cu = bdku, 



(7) 



where ku is the shear stiffness per unit area of the lamina and hence the effective shear stress r is 
defined: 



T = kuSu- 

The work done in the shearing process is given by the expression: 



(8) 



Us 



(5n(a) 



Fuia')d 



t 



cos j3 



[sin(a' ~ (3) + sin/?] 



kubdt 
2 

kubdt^ 



2 r 



sin(Q; — /?) + sin /3 



cos /? 



(9) 
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Figure 6: Two internal laminae of the laminated composite model, the shaded region shows the 
interlamina region which is exaggerated in scale for clarity. Dilation and shearing forces with their 
corresponding displacements arc given by: Fi and Fu with i5i and Su respectively. The highlighted 
section shows the lengths AB and BC that directly relate to 5i and 5u respectively. 




(a) sgn(<5c) = sgn(<5M) (b) sgn((Jc) = -sgn{(5M) 



Figure 7: Picccwisc linear force versus displacement model applied for interlamina shearing: (a) a 
linear-softening response which is more representative of a fracture model; (b) a linear-hardening 
response which is more appropriate for materials that show post-yield strength. 



for a ^ ac, or: 



2 
+ 

kubdt^ 
2cos2/3 



sin(a' — /?) — sin(aM — /3) 



sin(Q:c - P) - sin(aM - ,5) 



[sin(ac - /3) + sin/3]^ 



sin(ac — /S) + sin/3 
cos^ /3 

sin(ac — (3) + sin /3 



sin(Q:c - - sin(aM - /3) 



cos(a' -/?) da' I 

sin^(Q: — /?) 



(10) 



— sin^(ac — /?) + 2 sin(aM — P)[sm{ac — P) — sin(a — /?)] 



S{a), 



beyond the proportionality limit where a = [acoiul- However, if a > Qm and sgn(5c) = 
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(a) Mode I 



(b) Mode II 



(c) Mode III 



Figure 8: Fracture modes. Mode I is tearing; Mode II is shearing and Mode III is scissoring. In 
the current model, only Mode II is relevant. 



sgn((5M) > the shear resistance force vanishes and the expression for Us becomes: 
kubdt'^ 



Us 



C{ac) 

sin(Q;' — /3) — sin(aM 



sin(Q;c 



sin(Q;M 



sin(ac — /3) + sin^ 



cos^ /? 



cos(a' - /?) da' 



2 cos2 /3 
kiihdt^ 



sin/3 [sin/3 + sin(Q!c — /3) + sin(aM — /3)] + sin(ac — P) sin(aM — 13) 



(11) 



There would still be the potential for frictional forces to resist shear even though 5 > 5m and 
the interlamina region has lost all shear strength. However, this effect is currently neglected and 
hence the current model would tend to underestimate the true strength to some extent. 



3.1.3 Remaining energy contributions 

As in the geological model, the strain energy stored in bending can be taken from a pair of 
rotational springs of stiffness c: 

Ub = CO?. (12) 

The stiffness of the rotational springs is related differently from the geological model as the ex- 
pression for that model contained the overburden pressure q [221 . Since the bending energy should 
strictly relate to curvature k, where: 

.6/2 , 

Ub = 2 t:EIk^ dx, (13) 

J~b/2 2 



with X defining the domain of one bending corner and k as the rate of change of the kink band 
angle a over the kink band width b, as represented in Figure |9l thus: 

,. ,, 6/2 




(a) (b) (c) 



Figure 9: Bending of a lamina: (a) definition of x; (b) idealized case; (b) actual case. Curvature 
K is defined as the total angle change 2a over the effective length of the band 2b, hence kk, a/b. 
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««-^C«-. (14) 

Hence, the rotational stiffness c is related to the flexural rigidity EI of a lamina with E being its 
the Young's modulus in the axial direction and its second moment of area / = dt^ / 12. The strain 
energy per layer associated with the in-line spring of stiffness k is hence given by: 

Um = IkSl, (15) 

where Sa is the axial displacement of the springs. The in-line spring stiffness k — Edt/L for a 
single lamina with L being the length of the strut. The work done by the external load can be 
taken simply as the sum of the displacement of the in-line springs da and from the band deforming 
multiplied by the axial load P, which can be defined as the axial pressure p multiplied by the 
cross-sectional area of a lamina, dt: 

PA = pdt [5a + 6(1 - cos a)] . (16) 
3.1.4 Total potential energy functions 

The total potential energy V is given by the sum of the strain energies from bending Ub, the in-line 
springs Um, interlaminar dilation Ud and shearing Us, minus the work done PA, thus: 

V = Ub + U,n + UD + Us~ PA. (17) 

Since the dilation terms are assumed to be linearly elastic throughout their loading history, the 
total potential energy per axially loaded lamina takes three forms: 

1. The case where a = [0,ac], so F = V'", i.e. linearly elastic in shear. 

2. The case where a > ac, so V = V^, i.e. the secondary shear stiffness is either a smaller 
positive value than the primary shear stiffness or a negative value. 

3. The case where a > au and sgn((5c) — sgn((5M) > 0, so V = V^, i.e. no shear stiff'ness, 
which only occurs if the secondary shear stiffness is negative. 



These forms of the total potential energy are given by the expressions: 

V^ = y^ + ^£(a), V^^V^+'^Sia), = + '^S{au), (18) 
where is given by: 



yi - Ml + Edt^a^ ^ kibdt^ \^ cos(a - 



2 126 2 



1 - 



cos^ 



- pdt [6a + 6(1 - cos a)] . (19) 



The total potential energy functions are nondimensionalized by dividing through by kt^ and can 
be re-expressed in terms of rescaled parameters: 

V^ = V' + t^C{a), V^ = V'+'^S{a), = + ^Siau), (20) 



where: 



~_6 ... _pd _pL ~ _ Ed _ L ~ _ kidt _ kiL ~ _ kudt _ kuL 
~? ~Uk~m' "~~^~~B~' 



(21) 
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3.2 Equilibrium equations 

The equilibrium equations are defined by the condition of stationary potential energy with respect 
to the end-shortening Sa, the kink band angle a and the kink band width b; these can be written 
in nondimensional terms, thus: 



p = S, 

p = hla + hiJa 
p = kih + fell Jb 



2Dc 



IP' sin a 
Da' 
IP (1 — cos a) 



(22) 
(23) 

(24) 



Equation p2|) defines the pre-kinking fundamental equilibrium path that accounts for pure com- 
pression of the in-line springs of stiffness k. Equations (|23p -(|24 l) define the post-instability states 
for the non-trivial kink band deformations; equating them allows the kink band width h to be 
evaluated analytically: 

( ~ ^ 1/2 

~ I Da [2/ sina + a/ (1 — cosa)] I 



[ h (h - la) + hi (Jb - Ja 

The expressions for /q, and /f, are given in detail thus: 



(25) 



Ir, 



1 



cos(a — /3) 
cos /3 



sin(Q; — /3) 
sin a cos /? ' 



1 



2 (1 — cos a) 



cos(a — /3) 



cos /? 



(26) 



where these expressions apply for the entire range of a. However, the expressions for Ja and Jb 
change for each form of the total potential energy function; for V = V^, the expressions are: 



cos(q; — /?) [sin(Q; — /?) -I- sin/3] [sin(Q: — /3) -|- sin/3]^ 

^a — : 9~7^ ; ^b 



sin a cos^ /3 



2(1 — cos a) cos^ /3 



(27) 



for V = V^: 



cos(a — /3) 
sm a cos p 



sin(Q;c — /3) + sin/3 



Jb = 



sin(Q;c - 
sin(ac — P) + sin /3 



) - sin(Q;M - 
sin(Q;c — f3) + sin /? 



sin(a — /3) — sin(Q;M — P)] 



2(1 — cos a) cos^ (3 
^ sin^(a — /?) — sin^(ac — /3) + 2sin(aM — /3)[sin(ac — P) ~ sm{a — /3)] 
sin(Q;c — (3) — sin(Q;M — /3) 



(28) 



and for V ^V"^: 



0, Jb = 



sin/3 [sin /3 -|- sin(ac — /3) + sin(Q!M — /3)] -I- sin(ac — /3) sin(Q;M — /3) 
2 cos^ /3(1 — cos a) 



(29) 



The initial limiting case where a gives 6 — > cx) and p fcitan^/3 + fcn- The result for b 
suggesting that the kink band is initially prevalent throughout the structure and the result for p 
showing that the critical load depends primarily on the shear stiffness with a smaller contribution 
from the dilation stiffness that relates to /3. This reproduces similar results from the literature 
where the critical stress is related to the shear modulus [H |21 [13; it also refiects a significant 
difference from the geological model which has an infinite critical load and where the kink band 
width grows from zero length |22) . 
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Rod fibre diameter: t = 7 x 10^"^ mm 

Overall rod length L = 76 mm 

Longitudinal Young's modulus: E = En = 130.76 kN/mm^ 

Transverse Young's modulus: £'22 = 10.40 kN/mm^ 

Shear modulus (initial to final): G12 = 6.03 — > 0.68 kN/mm-^ over 4% shear strain 

Nondimensional stiffness quantities: 

Flexural rigidity: D = L/{\2t) 

Dilation stiffness: ki — E22L / {Ent) 

Shearing stiffness: fcn = Gi2L/{Eiit) 

Table 1: Properties used in the validation study to compare the current model with experiments 
presented in [S]. 

4 Numerical investigations 

4.1 Validation against experimental results 

Results from the current model are initially compared with published experiments on circular 
cylindrical composite rods with confined ends that exhibited kink bands under axial compression 
[5] . Although the model is formulated for flat rectangular laminae, the lamina thickness t can be 
perceived to be equivalent to the diameter of an individual fibre rod. This aids the comparison 
between the current model and the experiments such that both loading levels and the kink band 
width can be compared; a similar approach was employed in \V2\. 

The dimensions of the overall sample had a diameter of 8.255 mm with the relevant properties 
given in Table [TJ Note that the breadth d is not given since it cancels in all the relevant nondi- 
mensional quantities. The sample comprised ICI APC-2/AS4 composite fibres. Since the sample 
was cylindrical, the system in [51 was presented in terms of a cylindrical polar coordinate system 
with xi and X2 being the longitudinal and the radial coordinates respectively, as shown in Figure 
1101 For the tests presented to measure the change in shear modulus G12, there was no plateau 







; 8.255 mm ( ] 


( 


1 — ►.x, : 


1 \ 

1 X2 




1 

76 mm 1 



Figure 10: Representation of the experimental sample in [5]. The rod comprised ICI APC-2/AS4 
composite fibres with properties as given in Table [T] The sample was confined such that there was 
negligible transverse compression but also that global buckling was not an issue. 

shown in the test data, see Figures A3 and A5 in [S] . In the current study, it is therefore assumed 
that the piecewise linear model for the shear stiffness reflect the initial and flnal values found 
in the experiments; hence sgn((5M) — ^sgn((5c), i.e. a linear-hardening model is implemented as 
represented in Figure (T^b). 

The critical shear angle, 7c, which is effectively equal to the so-called engineering shear strain, 
for the piecewise linear idealization, is the angle beyond which the shear stiffness is replaced by 
a secondary smaller value; this is estimated from the aforementioned graphs in Figures A3 and 
A5 in to be 0.012 rad (or 0.69°). The shear angle can be expressed in terms of the kink band 
orientation angle /3 and the kink band angle, a, such that: 

^ii(ac) sin(ac -/3) +sin/3 
Si(ac)+t cos(ac-/3) 

Given that /3 is assumed to remain constant during deformation, the critical kink band angle ac 
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can be found by rearranging (j30[) . thus: 

tan 7c cos (ac — /3) — sin (ac — /?) — sin/3 = 0, (31) 

and solving for ac- This is achieved by substituting the critical shear angle 7c from above and 
the kink band orientation angle from [5], where /3 was reported to lie between 12° 16° to the X2 
direction; for the specified values, ac is approximately equal to 7c (ac = 0.0120 — > 0.0121 rad). 
To obtain the correct final shear modulus (see Table [T]) the value of Su/t ~ —0.094 is used 
such that the ratios between the initial and final values of the shear stiffness reflect the reported 
experimental data. 

Figure [Tl] shows numerical results from the current model using the properties defined in 




15 




Increasing /? 



80 ■ 
60 

40 

20 ■ 



0.2 



0.4 




0.01 



0.02 0.03 

Sa/t 



0.04 0.05 



Figure 11: Nondimensional plots of load p versus (a) total end-shortening A and (b) kink band 
angle a (rad) ; (c) kink band width b versus kink band angle a (rad). Range of /3 = 12° 16°. 
(d) Piecewise linear-hardening relationship of the effective shear stress r (N/mm^) versus the 
normalized shearing displacement Su/t for /3 = 16°. Properties of ICI APC-2/AS4 composite 
fibres and configuration and the range for /3 were taken from [5] . 

Table [T] with /3 values as found in the published results. Note that the nondimensional total 
end-shortening A is defined thus: 

A = ^ + 6(1 -cos a). (32) 

The actual kink band widths in the 5 tests were reported to range from 11 to 36 fibre diameters 
(directly corresponding to b in the current model) and the compressive strengths were found to 
average at 1.119 kN/mm^ with a standard deviation of 0.043 kN/mm^ (directly corresponding 
to p in the current model). The results from the current model show highly unstable snap-back 
and hence the critical load would never be reached realistically, see Figure [TTT a) and (b); a well 
established feature for systems of this type [4]. For comparison purposes, the pressure p is taken 
at the point at which the structure stabilizes and reaches a plateau; for the range of the f3 angles 
considered, the nondimensional stabilization pressure p ranges from 80.5 — ?> 82.3 which converts to 
an actual stabilization pressure p ranging from 0.969 kN/mm^ 0.992 kN/mm^: an error against 
the average from the experimental results of between 11% ^ 13%, which is sufficiently small to 
offer encouragement. 
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Kink band angle 



Case: l3 = 12° 



Case: /3 = 16° 



a = l3 
a ^213 



b= 10.4 



b 



17.0 



b= 10.2 



b= 19.5 



Table 2: Nondimensional kink band width values from the current model at different stages of 
deformation. The conditions a — jS and a — 2/3 are the points where the dilation within the band 
are effectively maximized and minimized respectively; experiments in [5j reported 6 = 11 — >■ 36. 

Of further interest is the comparison for the kink band width between the tests and the current 
model. Observing the graph shown in Figure [TTJc) , as the kink band angle a increases, initially 
the nondimensional kink band width b falls from a large value to a small value, approximately 
5.6 when a — 0.024 rad(w 1.4°). As a increases further, the kink band width begins to increase 
slowly; see Table [5] for details of some key points. According to the sequence described in Figures 
IHb)-(d) the kink band itself maximizes dilation when a = P, minimizes it when a = 2/3 and locks 
up when a > 2(3. The results of the current model, particularly when a = 2(3, lie at the lower end 
of the range of observed values of the band widths from the published experiments. This seems 
sensible given that the lock-up condition used, where a = 2/?, represents a lower bound [22], which 
implies that the current model would also tend to predict lower bound kink band widths. Hence, 
the results from the comparisons between the current model and the published experiments [5] are 
highly encouraging; they offer very good quantitative agreement for the loading and the geometric 
deformation - key quantities that define the kink band phenomenon. 

4.2 Parametric studies and discussion 

The favourable comparisons between the current model with the published experiments in [5] 
imply that the fundamental physics of the system are captured by the current approach. The 
study is therefore extended to present a series of model parametric variations to establish their 
relative effects. The basic geometric and material configuration is identical to that used in the 
validation study presented in Table [T] Material and geometric parameters are varied individually, 
while maintaining the remaining ones at their original values. The parameters that are varied 
are the kink band orientation angle (3, the critical kink band angle, ac, the composite direct and 
shear moduli, En, E22 and G12, and the shape of the piecewise linear relationship for shear. 

4.2.1 Orientation angle 

In the current model, the orientation angle (3 needs to be fixed a priori, hence the effects of different 
starting conditions for the model need to be established. Increasing (3 from 10° to 30°, a range 
that is representative of laminate experiments in the literature [5] [24] , leads to a stiffer response 
for increasing a, as shown in Figure [T2l with the pressure capacity for (3 — 30° being more than 
double the capacity for (3 — 10° for values oi a < (3. The graphs in Figure [13] raise an interesting 
point about the response particularly when (3 (^ 22.5°, which seems to define a boundary where 
the kink band width b loses its monotonically increasing property after it initially troughs for a 
small value of a, which was identified as approximately 1.4° in H4.1[ Both graphs show that the 
kink band width at the lower bound lock-up condition a = 2/3 temporarily peak when (3 = 22.5°. 
For higher orientation angles the kink band width b in fact peaks beyond a = 1.4°, then troughs 
and then resumes the monotonic rise as seen for /3 ^ 22.5°. Moreover, this also explains the 
reason why the stabilization pressure increases for (3 > 22.5°, as shown in Figure [TWa) . since the 
pressure has an inverse square relationship with the kink band width as shown in the equilibrium 
equations (|23]) and ((24]) . The graphs presented in Figure [14] attribute this loss of monotonicity in 
b (beyond a = 1.4°) to the dominating influence of the dilation terms for larger (3, particularly 
in the region of maximum dilation where a ~ (3. In the flrst instance, it should be recalled that 
when (3 is larger the potential maximum dilation displacement 61 is also larger relative to Su when 
a = (3. Figure [Mj^a) shows that the maximum of the dilation term ki{Ib — la) from the expression 
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Figure 13: (a) Nondimensional kink band widths b versus the kink band angle a (rad) for a range 
of orientation angles /3 = 10° ^ 30°; circles mark the the lower bound lock-up condition a = 2/3. 
(b) Values of nondimensional kink band widths b and applied axial pressure p at the lower bound 
lock-up condition. 



for &, i.e. equation (1251) . increases substantially with (3 whereas Figure ITW b) shows only very 
marginal changes in the respective shear term kii{Ji, — Ja). The numerator in equation (j25|) . &num, 
which represents the influence of bending, is independent of /3, as shown in Figure I14f c) but the 
respective denominator, 6dcn, shows that the dilation term influences the values significantly for 
the higher f3 values, as shown in Figure fMld). Once a gradually increases above /? the dilation 
displacement progressively reduces and the shear term begins to dominate with the result that the 
kink band width resumes growth and lock-up occurs. This effect is similar to that found in the 
geological model with the introduction of the foundation spring of stiffness kj [22], as shown in 
Figured] the kink band width was also found to plateau with higher foundation stiffnesses. It is 
worth noting that if destiffening in the constitutive law for dilation was introduced that the effect 
found in the present case would be generally less pronounced. 

4.2.2 Critical shear angle and modulus 

Increasing the critical kink band angle from ac = 0.69° — ?■ 0.96° (with a fixed limiting displacement 
Su/t = —0.094 as before) shows an increase in the critical shear stress before destiffening occurs - 
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Figure 14: Graphs of various terms from the expression for the nondimensional kink band width 
b, equation ([^5]) for /3 — 10° — ?> 30° versus the kink band angle a (rad). (a)-(b) Plots of dilation 
and shear terms respectively, (c)-(d) Plots of the numerator and denominator of the b expression. 



see Figure llSf e) - and leads to a monotonia increase of the axial pressure p and the minimum kink 
band width b - see Figuresirsla) and (c). A subtly different pattern is observed in FiguresfTSlb) and 
(d) where trends for increasing the initial shear modulus — 6.03kN/mm^ — > 8.45kN/mm'^), 

lead to higher stabilization pressures but smaller minimum kink band widths. These are logical 
results since the effect of increasing the critical kink band angle will lead to a later destabilization 
in shear and hence increase the load and band width; the increase in the initial shear modulus 
increases the resistance against shearing - the process of kink banding therefore requires more 
axial pressure to overcome this increased stiffness. However, the increased shear stiffness reduces 
the kink band width since there is a greater resistance to that type of deformation. 



4.2.3 Hardening and softening in shear 

The variation in the piecewise linear model for the shearing response is now discussed. The 
constitutive behaviour, Fu versus has been hitherto assumed to be a linear-hardening law 
which corresponded with the data from the literature used in the validation exercise. Figure 1161 
shows results for different secondary slopes while they remain positive (a linear-hardening law). 
Figure [17] shows results for reducing the secondary slope further such that they become negative 
(a linear-softening law). The results exhibit progressive behaviour; the reduced secondary slopes 
reduce the load carrying capacity but increase the kink band widths. This can be understood from 
the softening of the internal structure giving less resistance to the kink banding process, allowing 
for larger rotations and gross deformations. For the cases presented in Figure [T71 the negative 
secondary slope mimics the behaviour of a fracture process where the shear stiffness and strength 
has vanished and mode II fracture and crack propagation would occur. However, as described 
above, a similar pattern remains with the strength reducing and the band widths increasing for 
weaker properties in shear, which appears to be entirely logical. The detailed effects of crack 
propagation have been left for future work although recent work on buckling-driven delamination 
[25] has suggested that an analytical treatment of such effects may indeed be tractable. 
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Figure 15: Comparison of the response of the kink band formation for different ac values in (a), 
(c) and (e) for the range of ac — 0.69° ~^ 0.96° and initial G12 values from 6.03 kN/mm'^ — > 
8.45 kN/mm^ (in (b), (d) and (f). Nondimensional plots of (a)-(b) load p versus total end- 
shortening A, (c)-(d) kink band width b versus kink band angle a (rad). (e)-(f) show the rela- 
tionship of the effective shear stress r versus the normalized shearing displacement 611/ 1. Note 
that /3 = 12° throughout. 



4.2.4 Young's moduli 



Results for a two- fold increase in the axial modulus En suggest that this only has a marginal 
effect on the stabilization pressure (~ 1.5% increase), whereas increasing the lateral modulus £^22 
results in a significantly stiffer response; the system stabilizing to a smaller kink band width (see 
Figure ITSl) . These, again, are logical results since the effect of increasing the lateral modulus 
E22 increases the resistance against dilation; the process of kink banding therefore requiring more 
axial pressure to overcome this. Increasing the axial modulus increases the axial stiffness fc, which 
in turn effectively reduces the relative dilation and shear stiffnesses without affecting the relative 
bending stiffness, see the scaling relationships in equation (|21|) . Since bending is currently assumed 
to be purely linear, its relative effect becomes progressively more pronounced and then outweighs 
the reduced dilation and shear effects at large rotations. Obviously, if the bending was assumed 
to plateau due to plasticity, this effect would be limited. 
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Figure 16: Nondimensional plots of load p versus (a) total end-shortening A and (b) kink band 
angle a (rad); (c) kink band width b versus kink band angle a (rad); (d) Piecewise linear-hardening 
relationship of the effective shear stress r (N/mm^) versus the normalized shearing displacement 
5ii/t. Range of 5M/t ^ -0.435 -0.094 and P = 12°. 

5 Concluding remarks 

An analytical, nonlinear, potential energy based model for kink banding in compressed unidirec- 
tional laminated composite panels has been presented. Comparisons of results with published 
experiments from the literature suggest that very good agreement can be achieved from this rela- 
tively simple mechanical approach provided certain important characteristics are incorporated: 

1. Interlamina dilation and shearing: the kink band rotation naturally causes shearing and 
changes the gap between the laminae, the matrix within the composite needs to resist both 
of these displacements for the laminate to have integrity and significant structural strength. 

2. Bending energy: the resistance to rotation sets a length scale which in this case is the kink 
band width b. 

Linear constitutive relationships for the mechanisms of bending and dilation, and piecewise linear 
for the process of shearing, together with nonlinear geometric relationships have been applied. The 
approach has been successful such that the mechanical response captures the fundamental physics 
of kink banding and agrees with the experiments from [S^ in terms of kink band widths and loading 
levels without having to resort to sophisticated numerical or continuum formulations. Unlike the 
geological model [H] , where a relationship was derived for the band orientation (3 that was related 
to the overburden pressure q, in the current case the angle /3 has to be assumed a priori since, as 
far as the authors are aware, no satisfactory procedure for predicting /3 for composite laminates 
exists. For laminates, the magnitude of the orientation angle /3 has been largely attributed to 
the manufacturing process [H I15j . However, if the overburden pressure is considered to be the 
controlling parameter for the equivalent "manufacturing process" that keeps the geological layers 
behaving together, then future work on modelling the process of manufacturing composite lam- 
inates may bear fruit; an indication of the parameters that govern the orientation angle for the 
current case may be established. Although this is a shortcoming for the present model, the results 
from the parametric study are very encouraging with the trends appearing to be entirely logical. 

The current model can of course be used as a basis for further work. In the first instance, the 
formation of new kink bands can be investigated since work derived from similar approaches to 
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Figure 17: Nondimensional plots of load p versus (a) total end-shortening A and (b) kink band 
angle a (rad) ; (c) kink band width b versus kink band angle a (rad) ; (d) Piecewise linear-softening 
relationship of the effective shear stress r (N/mm^) versus the normalized shearing displacement 
Su/t. Range of Su/t = 0.096 ^ 0.482 and /3 12°. 

the current one exist for the geological model [23j[26]. In particular, in j23] the lock-up criterion 
was employed as the condition to introduce a new kink band; although it was assumed that the 
original kink band stops growing, the comparisons between theory and experiments were shown to 
be very good. The piecewise linear formulation applied currently for shear could also be extended 
to include dilation giving the possibility of mixed mode fracture [27\ for the first kink band. 
Moreover, loading cases that are more complex than uniform compression could be investigated; 
for example, numerical approaches have been developed in |17[ I28j with varying degrees of success 
to investigate the formation of kink bands where there is a combination of shear and compression. 
An additional complication in the combined loading case is that the kink band propagation across 
the sample tends to occur more gradually in contrast to the present case where the formation 
process is fast. 
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